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ABSTRACT 


A simple  relation  is  derived  between  the  thick  turret  yield  of  neutrons  produced 
following  ta,n)  reactions  in  compounds  or  homogeneous  mixtures  and  the  yields  in  the 
separate  constituents. 


It  is  necessary  to  assume  that  the  Bragg  law  of  additivity  of  stopping  powers 
holds  for  the  compound  und  to  make  the  approximation  that  the  stopping  pavers  of  ull 
materials  vary  with  energy  in  the  same  way  apart  from  eonstunt  factors. 

The  range  of  validity  of  the  second  hypothesis  is  examined  in  sane  detail.  It 
is  shown  to  be  far  superior  to  the  approximation  of  a stopping  paver  independent  of 
energy. 
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Introduction 


The  production  of  neutrons  from  actinide  a-particle  emitters  in  contact  with  light 
elements  is  a problem  often  encountered  in  present  nuclear  technology.  Usually  the  medium 
is  a complex  mixture  or  compound  and  even  ir  the  cross-sections  for  neutron  production  are 
known  as  a function  of  energy,  or  the  thick  target  yields  are  known  for  all  tile  constitu- 
ents, it  is  not  a straightforward  matter  to  evaluate  the  thick  target  yield  in  the 
compound  or  mixture.  It  is  the  purpose  of  this  article  to  point  out  one  practical  method 
of  doing  so. 

Method 

Consider  first  the  separate  constituents  having,  in  isolation,  n^  atoms/c.c.  and  the 
mixture  or  compound  having  n^  atoms  per  c.c.  of  each  of  the  J constituents.  The  total 
yields  of  neutrons  tier  a-partlcle  for  the  separate  mnterials  are 

Y - n t j OjlEldx  (t) 

o 

where  are  the  ranges  of  the  a-partic les  and  o^iE)  are  the  cross  sections  for  the  (a,n) 
reactions  at  energy  F..  This  can  be  written 


Oj(E)  < dx/dE T dE 


(2) 


where  E is  the  initial  energy  of  the  a-particle  and  (dx.  <IE)  is  the  reciprocal  of  the 
a J 

stopping  power,  more  usually  found  in  tables  us  dE/d(p ^x) , where  pj  is  the  density. 

The  yield  for  a mixture  or  compound  can  be  written 


mixture 


^ "j  / °j 

V J 


),E»  ,av'*'rt.u.re  « 


(S) 


The  problems  in  using  this  expression  arise  because  both  Oj(E)  anil  vary 

considerably  with  energy  and  it  is  also  not  apparent  that  values  of  Yj  alone  can  ever  be 

use<l  to  evaluate  Y . . if  OlEl  and  (dx/dE)  . are  not  both  known  as  a function  of 
mixture  j j 

energy.  One  method  of  using  the  cross  section  data  is  to  neglect  the  variation  of 
(dx/dE )m^xturp  with  energy  and  take  it  outside  the  integral  in  equation  (3).  This  is 
not  us  bad  an  approximation  os  at  first  sight  appears  because  0^(E)  is  zero  below  the 
threshold  of  endothermic  reactions.  An  approximation  of  much  wider  applicability  how- 
ever, which  doesn't  neglect  the  variation  of  dF./dx  with  energy  and  which  moreover 
enubles  the  thick  target  duta  to  be  used  alone,  is  us  follows. 

There  is  a wide  range  of  atomic  number  ami  energy  for  which  the  variation  of  dE/dx 
with  E is  similar  for  different  materials  and  one  can  write 


(dE/dxJj  = (dE/dx) L 


(4) 


where  is  s constant)  independent  of  energy , for  the  two  Materials  J and  1 where  1 
refers  to  some  standard  material. 

Alternatively  one  can  use  the  tabulated  values 


(dE/dpx)  j 


% 

Pj  *J1 


( dE/dpx) 4 


(5) 


where  p is  the  density. 

From  equation  (2),  using  equation  (4) 


Y 


J 


Oj(E)  (dJt/dEjj 


dE 


(6) 


To  get  any  further  it  is  necessary  to  invoke  the  Bragg  additivity  law  of  stopping  powers. 
Without  this  law  or  its  e<gii valent  nothing  further  can  be  done  even  with  perfect  know- 
ledge of  Oj(E),  Yj  and  (dE/dx)j  for  the  constituents. 

The  Bragg  additivity  law  can  be  written 
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Since 


^mixture 


N Yj  "j  AJ 

J 


and 


PJ 


= njAj/N 


where  Aj  is  the  atomic  weight  of  element  j and  N is  Avogadro's  nunfcer,  equation  (7)  can 
be  simplified  to 


( dE/dx) mixture  = Y (dE/dx)J  * "/"j 
J 


or  using  equation  (4) 


( dE/dx)  mixture  = <<*/<*>! 

J 
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Inserting  this  in  equation  (3) 


mixture 


_J !o_ 


Oj(E)  . (dx/dE)  jdE 


Enj  V"j 


(8) 


Referring  to  equation  (6)  the  term  under  the  Integral  is  seen  to  be  directly  related  to 
Yj  and  so 


mixture 


y v / yi 

L ",  J1  J / n 
J J j J 


Jl 


(9) 


Thus  we  have  expressed  Y^^^  as  a linear  combination  of  Y^'s  the  thick  target  yields 
of  the  constituents.  The  weighting  factor  Tor  each  Yj  is  n'  Kjj/V  and  insofar  as 
is  a constant,  independent  of  energy,  for  the  materials  J and  1 and  the  Bragg  additivity 
law  holds,  the  expression  (9)  is  exact.  No  knowledge  of  ctj(E)  is  required  if  the  thick 
target  yields  are  known. 

The  extent  to  which  It  j can  be  taken  as  a constant  will  now  be  examined  using 
published  tables  of  stopping  power  for  a-particles.  There  are  three  widely  used  compila- 
tions by  Williamson,  Boujot  and  Picard  1966,  Northcliffe  and  Schilling  1970  and  Ziegler 
1977.  The  trends  of  stopping  power  are  similar  in  all  three  compilations  but  they  differ 
in  detail.  The  tables  of  Williamson,  Boujot  and  Picard  have  been  used,  in  spite  of  being 
the  oldest,  because  the  data  are  presented  more  conveniently  for  our  purposes. 

We  have  evaluated 


dE 

dx 
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dE 
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dE 
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/ dE 

J,5MeV  dxOxygen,5MeV 


(10) 


That  is  the  stopping  powers  in  different  materials  at  different  energies  are  compared  with 
that  of  oxygen  at  5 MeV.  In  Figure  I are  plotted  contours  on  the  Energy  - Z plane 
showing  the  departure  of  P from  unity.  The  surface  P=1  corresponds  to  K.j  value  which  is 
constant  and  independent  of  energy  where  the  subscript.  1 in  this  case  refers  to  oxygen. 

The  contours  define  a saddle  shaped  surface  with  the  greatest  departures  from  unity 
occurring  at  lo*  energies.  It  is  apparent  that  there  is  a wide  area  over  which  approxi- 
mate constancy  of  Kjj  obtains.  For  instance  if  we  take  the  ilC#  contours  the  range  of 
constancy  of  within  these  limits  extends  for  Li  from  1.9  MeV  to  well  beyond  10  MeV 
whereas  for  Uranium  the  range  is  much  less,  7.1  MeV  to  3.4  MeV.  The  figures,  it  should 
be  remembered,  refer  to  the  range  of  variation  for  the  particular  elements  with  respect 
to  oxygen. 
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Figure  I may  also  be  used  to  decide  on  a best  energy  at  which  to  evaluate  for  a 
particular  case. 

Suppose  for  instance  we  have  6.2  MeV  a-particles  incident  on  Beryllium  oxide.  The 
Z =A  abscissa  in  Fig.  1 intersects  the  -1.0#  contour  at  6.2  MeV,  the  +2.5#  contour  at 
3.3  MeV,  the  contour  at  2.5  MeV  and  the  + 1Q#  contour  at  1.6  MeV.  For  the  best 
average  K between  6.2  MeV  and  3.3  MeV  the  energy  at  which  Z=4  intersects  the  contour 
(-l+2.5)/2  = 0.  7SpC  should  be  used  ( ie)  at  4.4  MeV.  Similarly  between  6.2  and  2.5  MeV 
the  point  where  Z=4  intersects  the  contour  (-l+5)/2  = 2#  should  be  used  (ie)  at  3.6  MeV. 
Finally  for  the  best  average  between  6.2  and  1.6  MeV  the  point  where  Z=4  intersects  the 
contour  (-l+10)/2  = 4.3#  should  be  used  ( ie)  at  2.7  MeV.  The  insensitivity  of  the  thick 
target  yields  to  the  low  energy  part  of  the  integral  in  equation  (6)  makes  it  unnecessary 
to  go  to  lower  energies.  The  extreme  limit  of  error  involved  in  choosing  the  three 
energy  ranges  above  is  (4.5-0.75)#  or  +1.9#.  When  it  is  remembered  that  for  Beryllium 
(dE/dpx)  varies  from  670  MeV/gnv/cm2  at  6.2  MeV  to  1680  MeV/gm/cm2  at  1.6  MeV  it  is 
apparent  that  the  use  of  constant  Kj^  values  as  compared  with  assuming  dE/dx  is  constant 
reduces  the  sensitivity  of  the  computed  yields  to  dE/dx  by  at  least  one  order  of  magni- 
tude. 

Equation  (9)  is  not  expressed  in  the  most  convenient  quantities  for  practical  use*. 
It  can  be  transformed  as  follows.  If,  instead  of  we  use  the  quantity  Cjj  where 


C 


jl 


— K 

PJ  J1 


(10) 


and  where  PA  and  pj  are  the  densities  of  materials  1 and  j.  From  equation  (5) 


(dE/dpx) j = Cjj  (dE/dpx) 1 


(II) 


Using  CJX  instead  of  in  equation  (9)  we  get 

Vlxtur.  ' 1 

J J 

If,  in  the  case  of  a compound,  the  molecule  contains  aj  atoms  of  atomic  weight  A ^ and 

there  are  n molecules/cc,  then 
m 

n'  = n a 

J m j 

Substituting  this  in  (12)  and  remembering  that  p^  = njAy^J  and  cancelling  nffl  and  N 
in  the  numerator  and  denominator. 


Y 


mixture 


J Aj  CJ1  YJ 


“j  AJ  CJ1 


(13) 


*1  am  indebted  to  Dr.  B.W.  Hooton  for  this  comment. 
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The  weighting  factor  for  Yj  is  non  a^  A^  C^  where  A^  is  just  proportional  to  the 
fraction  by  weight  of  element  J In  the  compound  or  mixture.  Cjj  Is  obtained  directly  from 
stopping  power  tables  using  equation  (II). 

In  conclusion  it  is  worth  aentioning  two  special  cases  of  relation  (13)  which  are 
often  encountered. 

(I)  A small  admixture  or  a foreign  material,  2,  is  present  in  an  otherwise  single 
component  material. 

From  ( 1 3) 


r , ai  Al  Y1  * *2  A2  c2.l  v2 

mixture  “ a A ♦ a,,  A C 

* 1 * * * f I 


if  therefore  a^Cj  ,/•»,*,  « I 


mixture 


Y * * * r Y 

) a)A)  L2,l  2 


(14) 


This  relation  can  be  deduced  directly  from  equation  (3)  without  Invoking  the  Bragg 
additivity  law  so  is  independent  of  it. 

(2)  If  one  component  of  a mixture  has  xero  neutron  yield  (ie)  Y.,  = 0 then 

Y - "I  A.  Y< 

mixture  a(  A(  ♦ a.,  A0  C,}  ( ' 

Cj  t is  still  Involved  in  spite  of  the  fact  that  Y.,  » 0. 
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